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Abstract 

The mean field type approach based on the self-consistent consideration of an 
effective field created by electron transfer is developed for a description of ther- 
modynamics of the Hubbard type models with an infinitely large on-site repulsion. 
This procedure, formulated according to 0], is an extension of the recently proposed 
generalized random phase approximation (GRPA). 

Within this approach, the thermodynamic properties of the two-sublattice pseudo- 
spin-electron model (the Hubbard model with local anharmonicity) are studied. 
Such a model can be used for a description of dielectric properties of YBaCuO-type 
superconductors along c-axis; pseudospins represent anharmonic motions of apical 
oxygens 04. It is shown that there are either phase transitions in the model with 
jumps of the mean values of a pseudospin and of electron concentration (in the 
u = const regime) or the phase separation (in the n = const regime). The phase 
transitions or phase separation are caused by pseudospin-pseudospin interaction as 
well as by electron transfer (the latter results in appearing of effective interaction be- 
tween pseudospins). The possibility of the ferroelectric type ordering of pseudospins 
is investigated. 

1 Introduction 

The Hubbard model @ is a fundamental one in the theory of strongly correlated electron 
systems. After the discovery of high-temperature superconductors it has been intensively 
studied, because the HTSC- compounds necessarily possess planes with strong correlation 
between electrons described by this model. Amongst the most laborious but promising 
approaches to study the Hubbard model there is a diagrammatic expansion with respect 
to atomic limit, which essentially differs from the diagrammatic technique for Fermi- 
systems ||. One of recently worked out developments is the GRPA proposed by Izyumov 
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et.al. in order to study the Hubbard model with Coulomb repulsion U = oo 0-0. 
Application of the GRPA approach to a pseudospin-electron system can be found in 
|J. However, this approach, being developed for investigation of correlation functions, 
is not a self-consistent theory. Therefore, lattice instabilities within it are treated as a 
result of singularities of susceptibilities that cannot suffice when the model undergoes the 
first order phase transition. The question how to calculate the thermodynamic quantities 
within the GRPA remains open as well. 

Returning to the possible approaches to a description of models with strong electron 
correlation, one should note that due to a success in obtaining exact results for simple 
models (Falicov-Kimball model [|10[ , , pseudospin-electron model with U = |L2[], etc), 
the approximation of the infinite dimension of space becomes very popular. However, 
there are some problems here, too. The single-site problem for the Hubbard model re- 
mains analytically unsolved; there are some difficulties with obtaining the higher order 
expansions in powers of 1/d. 

In this paper we use the mean field approximation, proposed in |l] for investigation 
of the t — J model. We show using as an example the Hubbard model that the diagram- 
matic series, summed over in this approximation expands the diagrammatic series of the 
GRPA, and, hence, we can talk about the self-consistent approach within the GRPA. The 
approximation is applied to the two-sublattice pseudospin-electron model which, being 



derived for description of ferroelectric type effects in HTSC of YBaCuO type (see [[13 
|TT|j), takes into account the two-sublattice structure of an apex oxygen subsystem in such 
superconductors. Coupling between the electron and pseudospin subsystems ( pseudospin 
variables are used to describe the local strongly anharmonic vibrations of apex oxygens) 



is performed via an interaction of the griiS* type proposed by Miiller [18 



The paper is composed as follows. In Section 2 the method of construction of mean 
field approximation within the GRPA scheme for the Hubbard model with U = oo is 
presented. Application of the method to the two-sublattice pseudospin-electron model 
is given in Section 3. Sections 4 and 5 contain the results of numeric calculations and 
conclusions, respectively. 



2 MFA approximation applied for Hubbard model 
with infinite Coulomb repulsion 

In this section we show how to obtain the MFA equations in example case 
of one-band Hubbard model and discuss the accuracy of that approach. The 

Hamiltonian of the Hubbard model without taking into account the influence of magnetic 
field reads 

i,s i i,j,s 

where /i is the chemical potential; U is a Coulomb repulsion of electrons on the same 
orbital; describes an electron transfer from the site j to the site « of a real lattice, 
is the annihilation operator of an electron with spin s (s takes two values - f and |), 
the electron number operator. 
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Usually near the atomic limit the model is considered using the basis of four states 
Ri) =\ n\,n\) 

| 1> =| 0, 0> 

(2) 



1) = 


0,0) 


2) = 


1,1) 


3) = 


0,1) 


4) = 


1,0) 



and Hubbard operators which are defined as follows X^ s =\ Ri)(Si \. 
Let us introduce the "bare" (zero-order) Matsubara Green's function as [[19 



(TXf(r)Xf(r')) Q _ £P9(r _ r /) f ±n±(s m ) , T > t' 



9ij[T T; " ll <[xf,xfj ± >o ^° \ ± n± (e m ) - 1, r<r', (3) 

where T means T-product, Xf 9 (r) = e TH °Xf q e~ TH ° is an operator in interaction 
representation and (. . )o is the average with the Gibbs distribution with the 
zero-order Hamiltonian. e p are the single-site energy levels of the zero-order 
Hamiltonian (at Uj = 0): 

e 1 = 0, e 2 = -2fM + U, e 3 = e 4 = -p. (4) 

The upper sign is chosen when the operators Xf 9 are of the Fermi- type (X 41 , X 23 , X 31 , 
X 24 and complex conjugate); otherwise, the lower sign should be taken. Correspondingly, 
[ ... ] + is an anticommutator, and [...]_ is a commutator. 

The results of perturbation theory and applying of Wick's thorem [|T9| , p0[] 
can be presented by elements which are denoted in the following way 

(xf) = ®, (xf x?r = arfey (XDo = C© _®) 

(x?x?xr n ) e = arfey^fe ( x r m )o = C®:®:®j 

9™(t - r')(Xr + Xf ) ee gff(r - r')Bf 9 = H&B* 

Let us write here the diagrammatic expansion for the average of X-operator 
in Heisenberg representation (TX/^r)) = (Te rH X 11 e~ TH ), i — 1: 



3 





Retaining the diagrams with only simple loops, we can neglect the diagrams which 
contain the Green's functions gfj(r — r'). In the limiting case U = oo, which we shall 
consider hereafter, we can neglect the Fermi-type Green's functions g?j, gfj and the Bose- 
type Green's function gfj as well. Discarding the mentioned class of diagrams, one can 
show that the results of the application of Wick's theorem do not depend on the so-called 
priority rules ||. 

Classifying the obtained result, we treat the diagrams 2, 4, 6, 7 as a "thickening" 
of the zero-order average (X u ) by the insets of the mean field type. The three last 
diagrams in [| should be discarded from the reasons presented above. There are also other 
two important type of diagrams: the diagrams 9, 10 renormalize the vertex parts of 
the Green's function of the initial diagram 3, whereas the diagrams like 11, 12, 13 can 
be interpreted as an complication of self-energy parts in the more simple diagram 3. 
After Fourier transformation all of the mentioned diagrams (unlike the diagrams 14-17) 
correspond to the one sum (or to a certain power of the one sum) over the wavevector k 
and frequency uj. As one can see from the figure, the diagrams 3, 5, 9-13 are those which 
are summed over within the GRPA when an equation of the Bete-Salpeter type is solved 

EM- 
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Let us introduce a set of diagrams corresponding to the generalized Hubbard-I ap- 
proximation (in the simplified version that is valid in the case of independent subbands 
||), denoting it by a "thickened" Green's function. 




(6) 



where = g pq (uj)B pq is the "dressed" Green's functions with the renormalized vertex 

part. We shall consider only the uniform phase in which all sites are equivalent. Then 
the equation for the vertex part B pq can be written as 





(7) 



Via pq we denoted the change of a spin projection in diagrams. Thus 41 = 31, and vice 
versa, 31 = 41. "Dressed" Green's functions with semiinvariants -C^9- are of g 41 (u;)-B 41 
or g 31 (uj)B 31 type. They creates loops with spin f or j, respectively. The series in the 
right hand side of the equation is analogous to the series which arises in the mean field 
theory for the Ising or Heisenberg model. 

Formally, expression (|7]) can be presented as 



m = ^Sp (B Pq e-^ff) - L pq , 



(8) 



where 



Hl ff = H + Q 41 Bt 1 + Q 31 Bf 1 =H +J2Q Pq B? q , {pg} = 41,31; 

{vq} 




(9) 



The renormalized Green's functions g pg (uj) is defined by the equation 



Hence 



9 Pq 



1 




1 



iu - e pq + QP q ion - e pq ' 



£ pq = £ p -e q 



Q 



pq 



(10) 
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That corresponds to the summation of diagrams of the 11, 12, 13 type. The expression 
( |TTD describes a renormalization of the single-electron excitation spectrum. Q pq plays the 
role of the irreducible self-energy part, by which the levels of the single-electron transitions 
are shifted. Using the renormalized Green's functions one can construct the loops Q pq 
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and L pq and, thereby, derive an expression for the irreducible part and equation for the 
auxiliary quantity L pq : 



Q 



l 

N 



7i+ (e) 



1 



Lm = ^E n + (e pq + t k B pq )-n + (e pq ) 



(12) 
(13) 



The expression for the free energy was obtained with making use of the analogy be- 
tween the diagrammatic series for B pq and series for the magnetization in the mean field 
theory for spin systems. Then so called one-tail diagrams correspond to simple loops in 
current approximation. One can show that the free energy can be written as 



F = /in — T In Sp (e~ 



■PH, 



eff 



Q pq B pq + J2D pq , 
{pg} {pq} 



D pq 



1 
2 



+ 




/w? ln(1+e " 



-P(&">+t k BPi) 



(14) 



One can test this result by introducing additional fields h pq B pq into Hamil- 
tonian ([I]). Then condition Jp- = B pq will be fulfilled, if one takes into 

PI \hpq= 

account (0) and (|I3|). 

In the paramagnetic phase B il = B 31 = B, Q A1 = Q 31 = Q. Hence, the full set of 
simultaneous equations defining the mean field approximation for thermodynamic quanti- 
ties of Hubbard model consists of the equations (§) and (|12) as well as of the condition 
for mean number of electrons per site 



n. 



(X 33 + X 44 ) 

One can show that standard condition for Hubbard operators 

(X 11 + X 33 + X 44 ) = (X 44 ) + B = 1 



(15) 



(16) 



is fulfilled automatically. 

Considered this result within the wide class of approximation for Hubbard 
model we noted that in the limit U oo it is equivalent to results obtained 

(see also []25|]) which describes the shift 



within the moment technique [23 



of Hubbard bands on value depended on spin a 

1 



B n 



nJl 



^ty(a£.a JO -(2n i _ (7 - 1)). 

#3 



(17) 



with n a = (n<j). One can show that there is clear correspondence between B a 
and Q pq : 

S41 _ „ D /S31 



Q 1 



-ntB 



Q 



-niB t . 
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But in MFA approach we have also equations for end parts for Green func- 
tions B pq (n a ) which is needed for consideration of processes where ordered 
phase takes place as well as the equation for free energy ([14]) required for 
investigation of phase transitions of first order. 



3 Two-sublattice Hubbard model with local anhar- 
monicity 

The Hamiltonian of the two-sublattice Hubbard model with local anharmonicity which 
includes also an interaction of electrons with local anharmonic vibrations, described within 
the pseudospin formalism, reads (|2l|) 



H = H e + H s + H e _ s + H s 



He = ~H E« + "a) + u EWi^ii + + E E ^ 



isa a jsa 



i 



i 

i,j a/3 



TT _ \ " \ "> jOl(3 qz qZ 

^s-s — 2 J ij J ia J j/3 



Here H e is the Hubbard Hamiltonian for two sublattices (subplanes in the case of the 
YBa2Cu307_,5 structure), H s + H s _ s is a pseudospin part of the Hamiltonian similar to 
the Mitsui model used in the ferroelectrics theory; if e _ s describes the interaction between 
the electron and pseudospin subsystems. The operators h s ia and S^ a act at the site i of 
the plane a (a=l,2). h describes the asymmetry of a two-minima potentials. 

Similarly to that performed for the standard Hubbard model (1X1) , let us introduce a 

T I f I T I 1" I 

cluster basis of states | Ri) =\ n-x, n^, n- 2 , n\ 2 ) ® \ S^, S- 2 ), or \ Ri) = \ n^, n^, n- 2 , n- 2 , Sf u SQ 
(here n\ a ,n\ a denote eigenvalues of the operators n\ a andn} a ), consisting of 64 vectors of 
state 

|i) =|o,o,o,o,t,T) 

(19) 

| 64) =| 1,1,1,1,|,|). 

Analysis of the thermodynamics of the model (O) without taking into account the electron 



transfer has been performed in |2T] , where the regions of ferroelectric pseudospin ordering 
were found and phase diagrams where built. Now we will investigate the influence of 
electron transfer on thermodynamic properties and phase transitions. Let us consider a 
pseudospin part of the Hamiltonian within the mean field approximation and an electron 
subsystem via the presented in the previous Section selfconsistent GRPA scheme. 

The cluster basis of states consists of a direct product of two single-sublattice | i,a,R) 
= | n ia , n| a , S? a ) (a = 1,2) ones. Each single-sublattice set of states includes eight 
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components at a given site i and sublattice a: 

1) =1 0,0,1) 

9\ —I 1 1 I \ 

(20) 





o,o, t) 


1 i) =1 


0,0,1) 


2)H 


1,1, T) 


|2)=| 


1,1,1) 


3)H 


o,i, T) 


|3)H 


0,1,1) 


4)=| 


i,o, T) 


|4)H 


1,0,1) 



Let us introduce the Hubbard operators X?f =\ i,a, R)(i,a, S | in this basis. Then 
the operators a^ a , Oiia, Sf a are related to the X-operators as 



— Y 14 _i_ \" 32 _i_ y 14 -I- \" 32 ^ — Y 13 — Y A2 -I- Y 13 _ :r 42 

4 



^^W-^f) (21) 

Then Hamiltonian (|T^) can be easily rewritten as 

H = H Q + H int + H s _ s , (22) 

where H = £ Ee p a X™ with 

p=l i,a 
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1,1 _ ,h 

1 T 2 , t-2 — ^2 
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2 2 2/i + C/ T |± 5 , £ f = -2/i + i7±| T5 (23) 



! — — ^ T V -T 2 ->- y, =2 "-i/»TU 

-3,3 _ 4,4 _ _,, _ ft , <? 3,3 _ 4,4 _ , ft fi 

-1 — fc l ~~ A 1 T 2 2 ' 2 — fc 2 ~~ A i - 1 -2~r2 



and 



#int = J! %(^ia + X ia + X ia + X ia) X { X ja + X jl + X jt + X fa) + 
ij,a 

+ X ~~ -^fa + X ia ~ X ia) X ( X ja ~ X jl + X ja ~ X ja)i ( 2 4) 

ija 

H s _ s is left unchanged. 

It should be noted that the model fll8| ) does not take into account the transfer of 
electrons between the planes. Therefore, if the interplane interaction between pseudospins 
is absent (Jl? = Jf^ = 0), and also owing to the system symmetry, the problem is reduced 
to two identical single-sublattice subproblems. 

Due to the relation 

(25) 

and also when a splitting of the mean field type Sf a Sjp = {S- a )Sjp + {Sj l3 )S- a — ^(S- a ){Sj l3 ) 
is used, the effective Hamiltonian of the problem reads 

H eff = ff +EE ~ E J««' (<S£>Sfa. (26) 

{pg} ic iaa' 

where J aa / = £ J£f*' . 

i 

When writing equations of self-consistent GRPA approach for the two-sublattice model 
( |18D in the approximation of independent subbands || and at U = oo, we can use ex- 
pressions of the previous Section. All the quantities Q pq , B pq , etc. get an auxiliary index 
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indicating the subplane a = 1,2. For instance, owing to the presence of pseudospin de- 
grees of freedom, in addition to B pq we also have B^. Then the summation means 

{pq} 

the summation over pq = 41, 41, 31, 31. An additional contribution to the free energy D pq 
we should change by D? 9 , where 

K q = ^ E ^ (1 + ^) - ^ E In (l + e'^^) ., (27) 

with 

fpq = F p _ F g _ rypq 

°a °a °a ^a • 

The expression for a free energy now reads 
F = nn - Tin Sp ( e -«*««>) - E + E ^ + \ E ^' TO ( 2 §) 



2 

{p?} {pq} aa' 



The system of equations to be solved 



hpq _ Sp(B p a q e- pH -ff) _ fpq 



Q p a q = ^Et k n + (eZ" + t k B^) 

k 



(29) 



4 

should be supplemented by a standard identity for Hubbard operators (J2 X^ p ) = 1> 

P =i 



by the equation for the pseudospin mean value (S^) = Sp(S^,e~^ eff )/Sp(e~P eff ), the 

2 

equation for the chemical potential (X 33 + + X 33 + X^)g ef . = n as well as by the 

a=l 

equalities (X 33 ) = (X 44 ), (X 33 ) = (X 44 ) in paramagnetic case. 

4 Numerical calculations and results 



In paper [£2|, where the pseudospin-electron model with zero transfer was studied, the 
existence of the phase separation resulting from the direct interaction between pseudospins 
was revealed. This effect takes place, as it was shown, in the certain region of the model 
parameters values in the regime of the fixed electron concentration (n = const). On the 
other hand, the simple pseudospin-electron model or the more complicated model flUf ) 
introduced for a description of possible phase instabilities in HTSCs of the YBaCuO type 
do not take into account explicitly the existence of CuO chains which are a reservoir 
supplying a charge carriers into the superconducting planes. 

Therefore, the model can be considered in two regimes: a) constant concentration, and 
b) constant chemical potential (the fixed value of the chemical potential /x is sustained by 
structural elements which are not taken into account within the model). 

The ground state diagram in absence of direct pseudospin-pseudospin interaction and 
transfer (J a p = 0,tij = 0) is shown in Fig.|I]. At T = the chemical potential coincide 
with the line ABCD at any value of concentration n. When temperature is different from 
zero, the chemical potential lies to the right of ABCD at n ~ and to the left of it at 
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|i,o,o,i,t,4.) D 

|i,o,i,o,t,;> 

10,1,0,1,14) 
10,1,1,0, t,l) 




|l,o,o,i4,t) 
|i,o,i,o,;,t) 

|0,l,0,l,|,t) 

|o,i, 1,0,1, t) 



|0,0,0,0,t,4) 




|0,0,0,0,|,t) 



A 



Figure 1: Ground state diagram in the fi-h plane in absence of direct pseudospin ( J a/3 = 0) 
interaction ( J a p = 0) and electron transfer (tij = 0). U = oo. The electron and pseudospin 
configurations are shown. 



n ~ 2. In the [i = const regime, the chemical potential is given by the line parallel to the 
ordinate axis. 

Analysis performed in the spirit of the mean field theory |21j has shown that in the 
regime of a fixed value of the chemical potential, the region of ferroelectric type ordering 
of pseudospins lies in the vicinity of the FCBE curve. This curve is a boundary between 
the ground states with the pseudospin configurations | |1) an d | IT}- If the interplane 
interaction J 12 is different from zero, the mean value rj = (S^ + S^) is the order parameter: 
r) = in the paraphase and rj ^ in the ferroelectric phase. At J 12 = we deal with 
two independent single-sublattice problems, and the order parameter is zero; however, a 
jump of the parameter £ = (Sf — S^) takes place. 

Figs-HH illustrate the influence of electron subsystem on a appearance of the phase 
diagrams. Let us note that in the Mitsui model (to which one can pass by putting g = 
and = 0) the phase diagram has a form of a curve symmetric with respect to h. At 
"switching on" the interaction between the electron and pseudospin subsystems g, the 
center of a region with ferroelectric ordering shifts from to g in accordance with the 
broken line FCBE position when the chemical potential is changed (Fig.|TJ) . Besides, the 
phase diagram becomes asymmetric with respect to the center of the pseudospin ordering 
region. As it has been already noted, at J12 = the ferroelectric phase does not exist 
(figs.||b,0b,|3|c), a discontinuous change of £ takes place. 

Let us now consider the influence of electron transfer on the phase transitions in 
the simplest approximation Hubbard-I. It has been shown in that at U —>■ 00 in a 
one-sublattice case for the Hubbard model with local anharmonicity, one has in Hubbard- 
I approximation two subbands, separated by the gap equal to g: e 31 (k) = £ 41 (k) = 
f +t k (X 44 + X n ) and e M (k) = e il (k) = -f + t k (X u + X n ) . In the considered in this 
paper approximation, these subbands are shifted by Q pq ; the mean values of the Hubbard 
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Figure 2: Dependence of temperature of the phase transition T c on the parameter h at 
different values of interaction parameters Jn and J 12 in the regime \i = const = 0. Thick 
lines correspond to the case Uj/g = 0.2; thin lines correspond to the case Uj = 0. a) 
Ju = J12 = g/%, b) Jn = g, J12 = 0. Phase transitions are of the second order (solid 
lines) or of the first order (dashed lines). 
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Figure 3: Dependence of temperature of the phase transition T c on the parameter h at 
different values of interaction parameters Jn and J± 2 in the regime ji = const = —g. 
Thick lines correspond to the case Uj/g = 0.2; thin lines correspond to the case Uj = 0. 
a) Jn = J12 = g/2, b) Jn = g, J i2 = 0. c) Jn = Ji 2 = 0,tij/g = 0.2. Phase transitions 
are of the second order (solid lines) or of the first order (dashed lines). 
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Figure 4: Dependence of the order parameter £ (a), band edges (b), thermodynamic 
potential (c), and electron concentration (d) on temperature in the regime /x = const = 
—g. The parameters values are: J\\ = J12 = 0, T/g = 0.1, t%jjg = 0.2, h/g = 0.01. 



operators are not calculated in the zero-order approximation but found self-consistently. 
In the case shown on Fig.||] (/i = 0), the chemical potential lies right in between of the 
mentioned bands. The number of charge carriers slightly changes with temperature and 
parameter h. The transfer can facilitate, but by itself is not able to cause the phase 
transition. 

At jj, = —g (Fig.^) the chemical potential lies near the edge of the lower band. Even 
unsignificant change in the spectrum can make the chemical potential leave the band 
(Fig.f|b) and provokes jump of parameter £ (Fig.|]a). In this case, interaction via electron 
transfer alone can cause the phase transition in the system. Fig.[|c indicates that the 
phase transition of the first order takes place. The behaviour of electron concentration 
here is quite specific (Fig.£|d), when it jumps from a certain value at T = up to when 
the chemical potential leaves the band. On the other hand, the existence of instable and 
metastable branches in Fig.|] indicates that the chemical potential is complicated function 
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Figure 5: Dependence of the order parameter r\ (a), band spectrum (b), and free energy 
(c) on concentration. The parameters values are: J\\ = J i2 = g/2, T/g = 0.1, Uj/g = 0.1, 
h/g = 0.5. 

of concentration. Therefore the situation when at a certain value of n the function n(n) 
has a jump becomes possible. In this case we have separation to two phases with different 
concentration, similarly to a seen in the P(V) diagram separation of vapour into gas and 
liquid in the gas-liquid phase transition. 

Fig|| shows that even though the solution of the system of equations indicates that 
the ferroelectric phase arose (Fig.|5]a), the separation takes place before that (Fig.|5|b). It is 
also indicated by a certain concavity in the free energy behaviour (tangent dashed lines in 
Fig.|5|c link the points with concentration values n 1; n 2 and n 3 , n 4 on which the separation 
take place). Hence, in the n = const regime we can talk only about the separation on the 
paraelectric and ferroelectric phases at concentration values n\ < n < n 2 and n 3 < n < 714. 
The ferroelectric phase exists (in the case considered in Fig||) in the concentration range 
n 2 < n < 713. Fig.||a and||b show the boundaries of the separated phase with and without 
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Figure 6: Regions of phase separation at different values of interaction J a/3 , T/g = 0.1. 
Thick lines correspond to the case tjj = 0; dashed lines correspond the the case ^ 0. 
a) ^ii = g, J12 = 0, t i:j /g = 0.2; b) J u = J 12 = g/2, tv/g = 0.1. 



transfer, positions of which are temperature dependent. At T = the region of phase 



separations reaches the sides of the 0</i/g<l,0< , n,*<2 [|22] rectangle. The transfer 
influence is ambiguous. In the first case (|6|a) it narrows the region of separations, whereas 
in the second case (|§b), the transfer only changes the shape of the region boundaries. It 
should be noted that in the presented figures, the high temperature is taken; therefore, 
at chosen values of model parameters and in the absence of direct pseudospin-pseudospin 
interaction, the separation does not exist. 



5 Conclusions 

The mean field approximation, being the simplest one in the theory of many-particle 
systems, is not trivial for the models with strong correlations of charge carriers. In this 
case, the effective field is created by transfer of electrons from site to site, not by regular 
averaging of a direct interaction. In contrast to the case of Fermi-systems, we have to 
solve here not only the equation for the self-energy part, but also some extra equations 
for the averages of the Hubbard operators |IJ]. In this paper we show that this approach 
means a selection of diagrams with only one sum over the wavevector k and frequency uj 
in expressions for the self-energy part and averages and is a generalization of the GRPA. 
Obtained shift of lower Hubbard band corresponds to results of the moment 
technique [23, 24, 25] in paraphase and in the limit U — > oo. The MFA equations 
for self-consistency parameters and expressions for thermodynamic functions allow one to 
analyze the equilibrium states and phase diagrams for the systems with strong correlation. 

In |l(J it has been shown how to obtain an equation of the dynamic mean field ap- 
proximation (d = oo) for the Falicov-Kimball model by means of summation of diagrams 
and selecting the single-site contributions. Nevertheless, the question of analytical con- 
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sideration of such contributions remains in general still open. Hence, we hope that a 
self-consistent approach, similar to the presented one, will give us a clue how to solve 
a single-site problem in the infinite dimensions for the another examples of models with 
short-range electron interaction. 

The mean field approximation applied to the two-sublattice model with local anhar- 
monicity indicates a possibility of the phase separation in the system. It has been shown 
in [^2| that the direct pseudospin-pseudospin interaction is sufficient to make a system di- 



vide into regions with different concentration of charge carriers. Instabilities with respect 
to phase separation were found in the ferromagnetic Condo model (see for example [26 ) . 
single-sublattice pseudospin-electron model with U = and without direct interaction 
between pseudospins f27j. From this point of view, our paper shows that the electron 
transfer not only does not destroy the separation regions, but facilitates process separa- 
tion and even is able to form separated phases by itself. This effect takes place when the 
chemical potential \i lies near the band edge, and the change of energy spectrum (due to 
ordering of pseudospins) in the /i = const regime moves the chemical potential into or out 
of the band. One can conclude that the electron (hole) transfer in the pseudospin-electron 
model creates an effective interaction between pseudospins. It is qualitatively different 
from the standard one JijS-Sj 0, but its role is similar. 

The problem of applicability of the approach proposed in this paper, as well as of 
the GRPA itself, requires a thorough analysis. The approach can be improved at the 
expense of complication of loop diagrams and vertices. Here the magnon Green function 
g 43 (k, uj) arise. Taking into account the magnon contributions is especially important at 
U 7^ oo, since without them the system of equations for the averages of the Hubbard 
operators becomes dependent on the priority rules. These problems are the subject of our 
subsequent studies. 
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